Mesoscopic mechanism of adiabatic charge transport 
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We consider adiabatic charge transport through mesoscopic metallic samples caused by a peri- 
odically changing external potential. We find that both the amplitude and the sign of the charge 
transferred through a sample per period are random sample specific quantities. The characteristic 
magnitude of the charge is determined by the quantum interference. 



X 



Let us apply an external potential 0(r, t), which is 
changing slowly and periodically in time to a metallic 
sample. This potential causes finite net charge Q trans- 
ported across the sample per period. This phenomenon 
known as adiabatic charge transport [1], has been inves- 
tigated in several papers[l-4] for closed systems at zero 
temperature. In the adiabatic approximation, quantum 
state of such a system is characterized by its ground state 
wave function corresponding to the instantaneous value 
of the external potential cj>(r, t). However, in real experi- 
mental situations exact eigenfunctions of electrons are ill 
denned: the electron energy levels are broadened due to 
inelastic processes at T ^ 0, and, in the case of an open 
system, are further broadened due to finite dwell time. 

In this article we present a theory of adiabatic charge 
transport in mesoscopic systems with " open geometries" . 
We demonstrate that at low T both the magnitude and 
the sign of Q are sample specific quantities. The typical 
value of Q in disordered (chaotic) systems turns out to be 
determined by quantum interference effects. We evaluate 
this value and find that it is much larger than the one in 
ballistic systems. This enhancement manifests of the well 
known fact, that at low temperatures all electronic char- 
acteristics of mesoscopic samples are extremely sensitive 
to changes in the scattering potential [5-8]. 

Let us start with a qualitative picture of the meso- 
scopic adiabatic charge transport. The wave functions 
of electrons in mesoscopic disordered systems are deter- 
mined by the particular realization of an impurity po- 
tential and exhibit sample specific spatial fluctuations. 
Therefore, the spatial electron density profile is chang- 
ing slowly in time, together with the external potential 
</>(r, t) . According to the continuity relation, variation of 
the charge density in time requires currents in the sys- 
tem. The question arises: What is the condition for a 
total charge transfer during one period to be nonzero? 
Let the pumping potential <fi(r, t) be characterized by a 
finite set of functions g(t) = {g a (t)}, a = 1, ..to, which 
are periodic with the same period to: 



(r,i) =<Kr,g(i)) =X>a(r) 5 a(i). 



(1) 



The time evolution of the set of functions g(i) represents 
a motion of a point in to dimensional space M.. Due 



to periodicity of </>(r, i), the trajectory C of this point is 
closed. The above mentioned currents lead to Q ^ 0, 
provided C encloses a finite area in M. This requires 
that to > 2. 

To calculate Q we will use the Keldysh technique for 
the Green function matrix equation [9] 

(idt -H - 0(r,i))G(r,r';i,i') - I e - ph ({G}) 

= 6(r-r')S(t-t'), (2) 

where G is a 2 x 2 matrix, G n = 0, G X2 = G A , G 2 \ = G R 
and G 2 2 = G K , with G R A K being the retarded, ad- 
vanced and Keldysh Green functions respectively. Ho is 
the Hamiltonian for electrons which includes impurity 
scattering potentials, and I e - P h denotes the electron- 
phonon collision integral. The solution of Eq. (2) can 
be expanded in term of the changing rate of the exter- 
nal field 4>(r,t), provided the time of the electron dif- 
fuses across the sample is shorter than the period of the 
external potential to. In general, G(r,r';t,t) depends 
on the value of the potential </>(r,t') at all the previous 
time t' < t. However, in the first order adiabatic ap- 
proximation G(r, r'; t, t) is determined only by the exter- 
nal potential and its first time derivative at time t, i.e. 
G(r,r';t,t) = &({<f>(r,t)},{<j>(r,t)}). The local time de- 
pendence in this approximation allows us to express the 
charge transfer Qi per period to along ith direction as 



dtTr{PiG K (t,t)} 



e I co l a {g)dg a (3) 



where Pi is ith component of the momentum operator, 
Li is the dimension of the sample along ith direction, 
i = x,y,z, Tr means integration over the space coordi- 
nates, and 



<4(g) 



— Tr{-^- 

dg a mL t 



G K (t,t)}. 



(4) 



Eq. (4) valids only in the leading order in fl = 2it/to- 
We have introduced differential 1 — form w^(g) on the 
m dimensional space Ai. Using Stoke's theorem one can 
convert the 1 — form integral along the trajectory C into 
the 2 — form integral over any surface S spanning C, 

Qi = e J dg a A dgpTr l a/3 (g), 
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<0(s) = {^-^--^-^-}Tr{-^-G K (t,t)}. (5) 
og a dgp dgp dg a 2mLj 

In an isolated quantum mechanical system, the 2-form 
7rig corresponds to the generalized adiabatic curvature 
of discrete eigenstates discussed in [3]. The wedge prod- 
uct is skew symmetric, i.e. dg a Adgp = —dgp Adg a . Note 
that Q in Eq.(5) doesn't depend on the spanning surface. 

In the zero order adiabatic approximation Keldysh 
Geen function G^(r, r', t) can still be expressed through 
retarded and advanced ones Gf ,j4 (r, r', t) and Fermi dis- 
tribution function n F {e) 



G 



A 



(G?-Gf)[l-2n F (e)] 



(6) 



where matrix elements of G e correspond to the instanta- 
neous hamiltonian Hq + <f>{v, g), 



G,[v.v'.l) = I di'G(r,r';t-^+J)e !et ' 



(7) 



and depend on time only through g(t). 

According to Eqs.(2), (7), the first order correction to 
the adiabatic approximation for G K can be written as 



5Gf(r,r',i)=i / dr"{[± - n F (e)](T RR - T 

dn F {e) 



de 



} (8) 



where for any p = (R, A), q — (R, A), 



T pq (r, r',r",t) = 2G p e (r, r", t)d t G^{v" , r', t)5 pq 

+Gl{v,v",t)d t <j){v",t)Gl{T"y,t) (9) 

The contribution of the first two terms in Eq.(8) to Qi 
can be neglected provided Li is much bigger than the 
elastic mean free path I, Substituting Eq.(8) into Eq.(5), 
we obtain 7T^(g), 



<s(g) 



d£i(P, g) dEi(a,g) 



dg a (t) dg p (t) 
^•(« ! g)=/^r 1 ^(r 1 )^(|--A) 



de 



*^Gf(r,r 1| t)Gf(r 1 y i t)| f 



(10) 



Due to the disorder, the charge Qi is a random sample- 
specific quantity. To characterize Qi we calculate its av- 
erage < Qi > and variance < (SQi) 2 >(<> stands for the 
averaging over realizations of the random potential). In 
the following we assume that L$~^> L z , L x , L y I, where 
Lff, is the electron dephasing length. In this case one can 
express < Q > and < (5Qi) 2 > through < 71"^ (g) > and 
variance D*^, (g, g') =< 6^(^)6^, p , (g') >. These 
quantities can be evaluated in a standard way Jl(| . Fol- 
lowing Eq. (10), < TT l a g > can be calculated in terms of 
diagrams in Fig. la, 




Fig. 1. a). Diagrams for < 7r^(g) >; b). Diagrams for 
IT^ ct , /3 ,(g, g'). Solid lines represent retarded(R) and ad- 
vanced^) electron Green functions, dashed lines denote 
the impurity averaging, wavy lines represent the exter- 
nal potential. Shaded triangles stand for T>(r, r'), shaded 
boxes for diffusons P(o>, r, r'; g, g'). In b)., the shaded 
triangles with two dashed lines represent Hikami boxes. 



< ^(g) >= J drdr>Mr)Mr') 



n dV(r, r') 



(11) 



where v is the mean density of states in the metal, \x F is 
the Fermi energy, Dq is the electron diffusion coefficient. 
Note that < n l a p > doesn't depend on the external field, 
i.e. on g in the leading order of (\>j [i F . 

^a(3a'f3' ^ s determined by diagrams in Fig. lb and can 
be presented in the form, 



dn F (e + w/2) dn F (e — uj/2) 



J dridridrdrV Q (ri)0 Q '(r' 1 )V ri P(ri, r) 



de de' 
d 2 A 

dgpdgp> 



(12) 



which explicitly demonstrates antisymmetry of the 
pumping perturbations. A is given as 

A = Re[D(u>, r, r'; g, g?)T>(w, r', r; g, g')]V rl 25(r' 1 , r') 

+V(lj, r, r'; g, g')V*(iu, r', r; g, g')V r; £>(ri, r). (13) 

2?(r,r') and T>(oj, r, r'; g, g') satisfy the following equa- 
tions 



A,V 2 D(r,r') = £(r-r') 



(14) 



{i[u + «ty(r)] + D V 2 }V(u;, r, r'; g, g') = S(r - r'), (15) 
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where <50(r,g,g') = 0(r, g') — <p(r,g'). For Eqs.(14), 
(15), we use usual boundary conditions: £>(r, r') = V(e — 
e',r, r') = when r or r' is at open boundaries [11, 12]. 

Let us now consider the sample sketched in Fig. 2 with 
two gates(labled by a — 1,2), biased with a.c. voltages 
of the same frequency and with a phase shift S = Si — S 2 , 



V a (t) = V x sm(flt + S a ). 



(16) 



In this case, m = 2. Let us assume that the potential 
induced in the metal by the voltages V a is screened with 
a screening length r much less than L x and 

g a (t) = sin(fii + 5 a ), 

Mr) = ^j?%o ( Z « *)*). (17) 

where C is the capacitance of the gate, W ^> ro is the 
width of the gate along z direction, are the z coor- 
dinates of the center of the gate 1,2; and 9{z) is the step 
function: 9(z > 0) = 1 while 9(z < 0) = 0. 




Fig. 2. Geometry of the sample. Shaded bars represent 
gates 1 and 2, crosses represent random scatters. 

To evaluate < Q z > in the leading order in Vo/mf <C 1, 
one has to substitute Eqs.(ll), (17) into Eqs.(5) and take 
into account that 

/ dg a Ndg Sj = [ dt{d t g a (t)gp(t) - d t gp(t)g a (t)}. 
Js Jo 

(18) 

Given the volume of the sample v = L x L y L z and the to- 
tal number of electrons inside the sample N we present 

< Q z > as 



< Qz >= f esm(S){ 



CV r 2 



) 2 N. 



(19) 



fo ~ 1 is a geometry dependent factor. 

To determine H % a p a ip> one has to solve Eq. (15). When 



CV rl D 



(20) 



this can be done by using perturbation theory with re- 
spect to g a {t). Keeping only the bilinear in g and 



g' contributions to T>(u, r, r'; g, g'), we express the so- 
lution of Eq.(15) in terms of T>° (u) , z , z') , the g inde- 
pendent solution of Eq. (15) with <j>(r, g) = 0. As 
usual, in the quasi-one dimensional case, we neglected 
x,y dependences of 2?°. Using Eq.(12), we then express 
n*0 a ,0, (g, g') in terms of 2?°. It is important to no- 
tice that n^a a ,a, is independent of g, g' in the leading 
order of CVoTq/vEt <C 1. Thus, according to Eq.(5), 
\/ < (3Qz) 2 > is proportional to the area S enclosed by 
the trajectory C. As a result for CVbrg/w -C Er, 



<{ 8Q z f>={es^5)\^fh{—). (21, 



T 



here the functions /i (77) has the following asymptotics, 



V < 1 

V > 1 



(22) 



In the limit CVoTq/v 3> Et 1 ^^(g) is a random quan- 
tity in 2-D space Ai of {g} with the " correlation length" 
(5g c | ~ E T v/CV r 2 (which is much less than unity). 
y/ < (SQ Z ) 2 > is determined by the amount of "flux" of 
random 71"^ (g) field that threads the loop C and should 
increase slower than the enclosed area S itself. This limit 
will be considered elsewhere. 

It follows from Eqs.(19), (21) that the standard devia- 
tion of Q z is much larger than its average, provided N is 
much bigger than G z 2 , (G z is the dimensionless conduc- 
tance in the z direction) and the amplitude of the external 
perturbation is not too large. In this case the value of Q z 
is entirely determined by quantum interference effects. In 
the opposite limit, Q z can be characterized by its average 
given by Eq.(19). 

According to Eq.(19), < Q z > is proportional to sin (5, 
i.e. this quantity changes sign together with 5 and van- 
ishes at 5 = 0. In fact this is as well correct for the charge 
transfer Q z in a specific sample of given realization of dis- 
orders, in the weak perturbation limit. One can see this 
from Eq.(5), taking into account that Tr l a/3 (g) is a random 
quantity independent of g when g is small. More gen- 
erally, the charge transfer Q changes sign when S — > — 5 
for arbitrary amplitude of the external perturbation al- 
though the simple 5 dependence in Eq. (19), (21) is valid 
only in the weak external field. Indeed, T-invariance re- 
quires that changing the direction of the trajectory C in 
the space M (which corresponds to changing sign of S 
for the case of two gates) should result in the change of 
sign of the charge: = —Q^ where <- J (^->-) corre- 

sponds to clockwise(countcrclockwise) motion along the 
same closed trajectory C. In the presence of magnetic 
field, this identity acquires a form 



Q^(H) = -Q^(-H). 



(23) 



Here (H) and (H) are charges which correspond 
to clockwise and counterclockwise motion along the same 
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closed trajectory C in the space M. and H is the external 
magnetic field. When the external potential is weak, one 
can neglect the g-dependence of field. In this limit, Q 
is proportional to the area S enclosed by trajectory C and 
TT l a p{g, ~ 0). Thus, Eq.(23) indicates that in the vicin- 
ity of g = T^apis) is an even function of H. At small 
amplitudes of oscillations of the external potential, this 
leads to Q < _ J (^_ > )(H) = H), provided without 

magnetic fields the system occupies a T-invariant state. 

As it is usual in mcsoscopic physics, the magnetic 
field dependence of Q z (H) exhibits random sample spe- 
cific fluctuations with a characteristic period AH ~ 
&o/L z L x . Here $o is the flux quanta (the magnetic 
field is applied along y direction). At high temperatures, 
when the dephasing length is shorter than the sample 
size, mesoscopic effects become exponentially small and 
Q z is determined by Eq.(19). 

Different mcsoscopic mechanism of the adiabatic 
charge transport has been discussed in[13]. Inelastic 
electron-phonon processes cause shifts of the center of 
mass of the electron wave functions and thus contribute 
to Qi. However, in the case of open geometry sam- 
ples, this contribution is small compared with Eq.(20) 
as (r e _p/ji?T) _1 <C 1, where r e - p h is the electron-phonon 
inelastic scattering time. 

In conclusion, we would like to stress that the d.c. cur- 
rent discussed above is proportional to the frequency of 
the oscillations of the external potential f2. This distin- 
guishes the considered above effect from the usual photo- 
voltaic effect. The latter effect in the low frequency limit 
is dominated by the relaxations of nonequilibrium distri- 
bution of electrons in the presence of external field via 
electron-phonon inelastic processes [13]. Such a mecha- 
nism leads to a randomly directed d.c. current propor- 
tional to the absorption rate of the external field even in 
the presence of a single pumping gate voltage. It means 
that the photovoltaic d.c. current is proportional to f2 2 , 
as the frequency dependence of the absorption rate of 
the external field is. In an open sample in the low fre- 
quency limit, the photovoltaic current is smaller than the 
adiabatic current by a factor Q/Et <C 1. 

It should be emphasized that the value of Q in a finite 
open mesoscopic system is not quantized. An approxi- 
mate quantization of Q can be achieved in a Coulomb 
blockade regime, for a pumping of the charge through a 
quantum dot, which is weakly connected with the source 
and drain[14,15,16]. Q turns out to be rather well quan- 
tized provided the dimensionless conductance of the de- 



vice is small. However, under these conditions, pumping 
is an entirely classical effect. It is not related to the quan- 
tum interference mechanism discussed in this paper. 
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